Orbital-selective conductance of Co adatom on the Pt(lll) surface 
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We propose an orbital-selective model for the transport and magnetic properties of the individual 
Co impurity deposited on the Pt(lll). Using the combination of the Anderson- type Hamiltonian 
and the Kubo's linear response theory we show that the magnetization and dl/dV spectrum of 
Co adatom are originated from the 3d states of the different symmetry. A textbook expression for 
the spin-dependent differential conductance provides a natural connection between magnetic and 
transport properties of Co/Pt(lll). We found that it is possible to detect and to manipulate the 
different 3d states of the Co adatom by tuning the tip-impurity distance in STM experiments. 

PACS numbers: 71.27.+a, 73.20.At, 75.20.Hr 



Spin-polarized scanning tunneling microscopy (SP- 
STM) and spectroscopy (SP-STS) are important tools 
of the modern physics for studying electronic and mag- 
netic properties of a transition metal adatom deposited 
on a surface [l| . Depending on the surface properties that 
can be insulating [2| or metallic |3| , magnetic [4| or non- 
magnetic Q , the STM conductance spectra demonstrate 
different types of the spin- dependent excitations such as 
Kondo resonance step- type excitations, quasiparti- 
cle peak The resulting dl/dV conductance spectra 
can be well reproduced in the framework of the local- 
ized (Heisenberg-type) [8|;12| or itinerant (Stoner-type) 
spin- model approaches [13 j. 

One of the fascinating results of the STM experiments 
is the possibility to manipulate and to control the indi- 
vidual spin as well as the simplest nanostructures. For 
example, the authors of Ref. 14j were able to construct 
magnetic chains of different lengths which displayed evi- 
dence of coupled-spin behavior. In turn Seratte et al. [HI 
have demonstrated a direct access to the spin direction 
of individual atoms. 

Recently, the 3d-orbital degree of freedom has at- 
tracted a lot of attention due to its important role in the 
massive variations of the magnetic and transport proper- 
ties of strongly correlated systems 2 12 1 ■ For instance, 
using X-ray reflectivity data the authors of Ref. [l8| have 
proposed a method to derive the orbital polarization of 
transition-metal oxide interfaces. In case of a transition 
metal adatom deposited on a surface the STM technique 
can be applied to probe and to manipulate the individ- 
ual 3d states. For this purpose one should use a system 
that demonstrates a strong variation of the electronic and 
magnetic properties depending on the symmetry of the 
specific 3d orbital. Such a variation results in a strong 
orbital polarization of the STM spectra. 

In this paper we propose the Co impurity deposited 
on the Pt(lll) surface to be one of the best candidate 
to probe the individual 3d states using STM technique 



FIG. 1. (Color online) Schematic representation of the com- 
puter simulation setup used in this work. The electrons hop 
from the tip to the surface via the impurity with the full 
Coulomb interaction matrix. 



since it demonstrates a strong orbital polarization of the 
dl/dV spectrum. Using the exact diagonalization of the 
many-body electronic Hamiltonian and the Kubo's lin- 
ear response theory we show that the main contributions 
to the magnetization curve and the peak of the spin- 
polarized dl/dV spectrum are due to the 3c? states of the 
different symmetry. Based on the obtained results we 
conclude that one can detect the different 3d states by 
varying the tip-impurity distance. 



Computer simulation setup. To study the transport 
properties of the Co/Pt(lll) system we construct the 
Anderson-type Hamiltonian that contains tip, impurity 
and surface states (Fig.l) 
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Here e^, and e p are energies of the impurity, tip 
and surface states, df a and &+. (c+ CT ) are the creation 
operators for impurity and tip (surface) electrons, V% q 
(Vi P ) is a hopping between impurity and tip (surface) 
states, Uijki is the Coulomb matrix element and the im- 
purity orbital index i (J, k, I) runs over the 3d states 
{xy, yz, 3z 2 — r 2 , xz, x 2 — y 2 ). In what follows we assume 
Vi q and Vi P are real. Using finite-temperature exact di- 
agonalization method [li} we calculated one-particle and 
two-particle correlation functions for 40 excited states, 
which is enough to work at the experimental tempera- 
ture of 0.3 K with the Boltzmann factor smaller than 
10~ 6 . 

The energies of the impurity and surface orbitals as 
well as the impurity-surface hybridization parameters can 
be estimated from the first-principles calculations by fit- 
ting procedure described in RefQ- The on-site Coulomb 
and intra-atomic exchange interactions were estimated 
to be U = 6.75 eV and J =0.9 eV 0. Since the spin- 
polarization of the tip effective orbital and its overlapping 
with the impurity wave functions are unknown, they are 
free parameters in our model that should be carefully 
chosen by using scientifically sounded arguments. 

The tip-impurity hopping process can be simulated 
by three z-oriented adatom orbitals that overlap with 
the tip orbitals. To take into account all the excita- 
tion paths between the tip and adatom we define one 
complex Wannier function centered on the lowermost tip 
atom. This Wannier orbital overlaps with yz, xz and 
3z 2 — r 2 orbitals of the cobalt impurity. We used a Slater- 
parametrization [201 ] to define the corresponding hopping 
integrals, V q ^ z 2_ r 2 = 2V qyz — 2V qxz . Thus the number 
of the adjustable parameters can be reduced to two items, 
the spin polarization of the tip and the tip-impurity hy- 
bridization. It is useful to define some bounds for these 
parameters. 

Since in the STS experiments Q the tip with out-of- 
plane polarization was used then the magnetization is an 
intrinsic property of the tip and it should be consider- 
ably larger than the external magnetic field ranging from 
-7 T to 7 T. The tip- impurity hybridization parameters 
strongly depend on the vertical position of the tip. It 
is natural to assume that the hybridization between tip 
and adatom should be considerably smaller than that be- 
tween adatom and surface orbitals (0.3 - 0.7 eV). Thus 
we have solved the model Hamiltonian with the tip- 




Magnetic field (T) 

FIG. 2. (Color online) Magnetization (in (j,b) of the impu- 
rity states as a function of the external magnetic field, B z 
calculated with A = 10 meV and V q3z 2_ r 2 = 100 meV. Mt 
is the magnetization of the tip. 

impurity hopping integral {V q ^ z 2_ r 2) ranging from 50 
meV (weak hybridization) to 150 meV (strong hybridiza- 
tion). Within the Slater-parametrization [2l| it corre- 
sponds to the tip-impurity distances of 4.1 A and 3.2 A. 
The spin polarization of the tip A = e q ^ — e g j, ranges 
from 10 meV (weak polarization) to 100 meV (strong 
polarization). In what follows we assume that the spin 
polarization of the tip is along z-axis. 

Magnetization of impurity states. According to the 
STM experiments [3j, the differential conductance spec- 
trum of the individual Co impurity as a function of the 
external magnetic field demonstrates an S-shaped curve 
form. Since there are no signs of hysteresis, the applied 
magnetic fields are not sufficient to induce magnetism. 
To reproduce this situation we have performed model cal- 
culations with the magnetic field B ranging from -25 T 
to 25 T and different values of V„Zz 2 -r' 1 ancl A- In all the 
cases we observed the similar behavior of the magnetiza- 
tion curves of the different orbitals. The typical example 
is presented in Fig. 2. One can see that all the 3d states 
can be classified with respect to the sensitivity to the 
magnetic filed. The orbitals of yz and xz symmetry are 
almost paramagnetic and show a weak response to the 
field. At the same time the slope of the 3z 2 — r 2 curve 
is the largest one, which is a result of a strong magnetic 
moment localization [7[. 

Spin-dependent conductance. Since the Hamiltonian 
Eq.(l) is diagonalized the obtained eigenvalues and eigen- 
functions can be used in order to calculate complex two- 
particle correlation functions. The most important one 
is the differential conductance that in the framework of 
the Kubo's linear response theory can be written 

dv Zui lo + L n — kj n t 

nn' 

where Z is the partition function, E n is the eigenvalue 
of the Hamiltonian Eq.(l) and j nn ' is a matrix element 
of the current operator. To avoid the extrapolation from 
Matsubara frequencies to real frequencies we use energy 
with a small imaginary part u> + iS. 

In our finite-cluster approach the current operator can 
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be derived from the quantum Hamiltonian, Eq.(l) by us- 
ing an approach proposed in Ref . }22| 
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The current operator has a symmetrized form and con- 
tains the average of the tip-impurity and impurity- 
surface contributions. Such an operator can be used 
to calculate the expectation value of the current pass- 
ing through the impurity in different regimes, i.e. in the 
linear response or in the steady state. A similar sym- 
metrization of the current contributions was used in the 
previous investigations, for instance in Ref. [23]] . 

The resulting differential conductance, Eq.(2) con- 
tains the sum over the spin and effective orbital in- 
dices. Let us first discuss the orbital dependence of 
the dl/dV spectrum. All the excitations of the same 
spin, inn'in'n can b e classified into the direct and indi- 
rect processes. The examples of the direct excitations 
are presented in Fig. 3(a). There are a tip-impurity- 

di,jb q <j\n)i a tip- 



surface excitation (process I) (n\c pa u, la u, ity w qa 
impurity-tip excitation (process II) (n|& + J <ii (T <i^ T &,j Cr |n) 
and a surface-impurity-surface excitation (process III) 
{n\Cp a di,ydf a Cp a \n) . Here q and p denote tip and surface 
effective orbitals, respectively. These direct excitations 
can be reduced to the form that contains the spin opera- 
tor of the corresponding impurity orbital, therefore, one 
should expect that these contributions are very sensitive 
to the applied magnetic field. 

It is not the case for the indirect excitations [Fig. 3(b)] 
described by 4-point correlation functions which are 
{n\Cp (T dj cr df (J b qi j\n) for the tip-impurity-impurity-surface 
transition (process IV), (n|6 + T dj (T (f^ T & 90 .|n) for the 
tip- impurity-impurity-tip transition (process V), and 
{n-\ c p>crd'jadf a Cp a \n) for the surface-impurity-impurity- 
surface excitation (process VI). Such processes are due to 
the full Coulomb interaction matrix taken into account 
for the impurity states. We believe that they play an im- 
portant role for simulating surface-assisted Hall effects 
[24] and left their analysis for a future investigation. 

Having analyzed the different types of excitations we 
are now in a position to connect the finite-cluster con- 
ductance with that experimentally observed [3j. For 
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FIG. 3. (Color online) Schematic representation of the di- 
rect (process I, II, III) and indirect (process IV, V and VI) 
excitations that contribute to the differential conductance. 



this purpose we should choose the terms that correspond 
to the tip-impurity-surface excitations (the process I for 
3z 2 — r 2 , yz and xz orbitals). However, these excitation 
processes are described by the two-particle Green's func- 
tion that is non-diagonal with respect to the site indices. 
It means that the spectral function can change the sign 
and, as a consequence, the resulting conductance can be 
negative. To overcome this problem we symmetrized the 
Green's function of the process I [5^]. Thus the spin- 
diagonal part of the STM conductance takes the form 
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As for the non-diagonal spin contributions to the dif- 
ferential conductance, Jnn'jn'n ( CT ^ a ') tnat can b e ca ^" 
culated by using Eq.(2) they correspond to the spin- flip- 
assisted tunneling processes (llj . Their contribution to 
the STM signal will be discussed below. 

The spin dependence of the conductance provides 
an access to the magnetic properties of the surface 
nanosystem. To describe the connection between conduc- 
tance and magnetic moment of the impurity a classical 
Heisenberg-type expression was proposed by Wortmann 
et al. 26]. Within this approach the spin-dependent 
part of the conductance is proportional to the product 
dlsTM/dV ~ MtMc'o, where Mr and Mco are the mag- 
netization of the tip and Co adatom, respectively. 



In turn, Fransson et al. [lfj and Delgado et al. [27 1 
have shown that the tunneling conductance contains a 
background contribution that is independent of the local 
spin, a dynamical conductance which is proportional to 
the local impurity spin, and a susceptibility contribution 
that corresponds to the spin-spin correlation function of 
the local spins. Using such an approach one can extract 
the local spin moment and analyze the spin excitations 
recorded in conductance measurements. 

We would like to stress that the previous considera- 
tions of the spin-dependent conductance are based on 
some approximations, which strongly limits their appli- 
cability. The crucial one is the absence of the quantum 
fluctuations between the tip, impurity and surface states. 
Such a mean-field approximation is probably valid for the 
system with localized magnetic moments. However, in 
case of the metallic surface the magnetic state of the im- 
purity can be considerably suppressed by the conduction 
electrons of the surface. Moreover, the localization of the 
magnetic moments can depend on the symmetry of the 
3d state and hybridization with the surface states. This is 
the case for the Co/Pt(lll) system where the yz and xz 
states are itinerant and 3z 2 — r 2 is localized 0. Thus the 
magnetic Hamiltonian and the conductance should con- 
tain the Heisenberg-type terms for the 3z 2 — r 2 states and 
Stoner-type terms for the yz (xz) states. These quantum 
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fluctuations and hybridization effects are naturally taken 
into account in our Anderson model approach. 

To study the interplay between transport and magnetic 
properties we consider the spin-matrix structure of the 
differential conductance, jjy = J2a<y' ~W L - ^ n the s P m 



space the tensor dl/dV can be expressed in the following 
form 



dl 



dl* 



(u, A, B)a, (5) 



dV y ' ' ' dV y ' ' ' dV 

where 1 is the unit spin matrix, a are the Pauli ma- 
trices, f£ = \Tr a §r and ^ = \Tr a %B are spin- 
averaged and spin-dependent differential conductivities. 
Varying A and B parameters one can simulate the dif- 
ferent tip-impurity magnetic configurations, for instance 
purely paramagnetic state (A = and B = 0), impurity 
paramagnetic state (A ^ and 5 = 0), ferromagnetic, 
antiferromagnetic or non-collinear states. We consider 
the proposed expression to be preferable in comparison 
with others since it preserves the spin polarization of the 
STM signal even when the magnetization of the adatom 
is small or zero. This is the case for Co/Pt(lll) where 
M Co ~0.05/z B (Fig. 2). 

In case of the collinear tip-impurity magnetic config- 
uration, B — (0,0, B z ) and in the atomic limit (U >> 
Vi q , Vi p ) the -j^r term can be associated with the magne- 
tization of the individual 3d states. The a;-component of 
the differential conductance, provides the informa- 
tion about transitions between different spin states of the 
adatom [27} . In contrast to the previous works we use the 
electronic Hamiltonian and, therefore, varying the on-site 
Coulomb interaction and/or hopping integrals it is pos- 
sible to study the spin dynamics of the system in the 
different regimes, for instance in itinerant, localized or 
mixed- valence states. 

We now turn to perform a quantitative analysis of the 
dl/dV spectrum of the Co/Pt(lll) system. The compo- 
nents of the spin-dependent conductance calculated with 
different values of the tip-impurity hybridization param- 
eters are shown in Fig. 4. One can see that in the case 
of weak hybridization [Fig. 4(a)] the spin- independent 
and longitudinal conductivities near the Fermi level are 
mainly due to the itinerant yz and xz states. The result- 
ing spectrum demonstrates a peak at the Fermi level and 
we believe that the model conductance calculated with 
these parameters corresponds to the experimental dl /dV 
spectrum measured by Meier et al. Q. Having calculated 
the conductance at temperatures ranging from K to 50 
K we did not find substantial changes of the spectrum. 
It is due to the fact that the width of the peak is much 
larger than the energy of the temperature excitations. 

In the opposite limit of the strong hybridization (Fig. 
4b) the localized 3z 2 — r 2 states give a considerable con- 
tribution to the resulting conductance and the system 
is in a mixed itinerant-localized state. The intensity of 
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FIG. 4. 



spin-dependent (^r- and d , 



(Color online) The spin-independent (-^S-) and 

dZf„ , dlf„ „ 

conductance components cal- 
culated in different regimes with B z — 8.6 T. (a) V q3z 2_ r 2 = 
50 meV and A = 10 meV. (b) V q3z 2_ r 2 = 150 meV and A = 
10 meV. The blue dashed and red dotted lines correspond to 
yz (xz) and 3z 2 — r 2 orbitals, respectively. The grey line is 
the experimental dlsTM /dV spectrum taken from Ref. 0]. 



the conductance at the Fermi level is considerably sup- 
pressed and there is a nearly insulating behavior of the 
dlo I dV spectrum that strongly differs from the weak hy- 
bridization case [Fig. 4(a)]. Thus we arrive at an impor- 
tant conclusion that it is possible to detect the different 
3d states of the Co impurity in the STM experiments 
by varying the tip- impurity distance between 4.1 A and 
3.2 A. For this purpose the STM technique should be 
combined with realistic many-body computer simulation 
methods. 

It is also interesting to analyze the transverse com- 
ponent of the conductance that is originated from the 
spin-flip excitation, for instance, {n\c^di^df,b q i\n) . The 
fact that using a Schrieffer- Wolff transformation the elec- 
tronic Hamiltonian Eq.(l) can be reduced to an effective 
tunneling spin Hamiltonian [^, 0, [ll| containing the spin- 
flip processes confirms that the original Hamiltonian also 
provides such excitations. The calculated yz (xz) and 
3z 2 — r 2 transverse contributions demonstrate remarkably 
different behavior. While the -^r spectrum of the itin- 
erant states is negative, the localized curve changes the 
sign near the Fermi level. The conductivities calculated 
in weak and in strong hybridization regimes are mainly 



associated with the yz and 3z 



states, respectively. 



The intensity of the spin-flip excitations varies strongly 
depending on the tip- impurity distance. At V q ^, z 2_ r 2 = 
150 meV such excitations strongly contribute to the re- 
sulting conductance since the system becomes more lo- 
calized and insulating. 

In conclusion, we numerically demonstrate that the 
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differential conductance of the Co/Pt(lll) system mea- 
sured in the STM experiments is orbital-selective. Vary- 
ing tip polarization and tip- impurity distance one can get 
the access to the localized and itinerant states of the Co 
impurity. The current SP-STS experiments operate with 
itinerant yz and xz states demonstrating a weak response 
to the applied magnetic field. The obtained results are 
useful for future scanning tunneling experiments aiming 
to manipulate individual 3d states of a transition metal 
adatom deposited on a surface. 
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